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Abstract 

We calculate corrections of order a {Za)Ep to hyperfine splitting in muonium generated by 
the gauge invariant set of diagrams with polarization insertions in the light by light scattering 
diagrams. This nonrecoil contribution turns out to be —2.63 Hz. The total contribution of all 
known corrections of order a^{Za.)Ep is equal to —4.28 Hz. 
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I. INTRODUCTION 



The hyperfine splitting in muonium is one of the intervals best studied both experimen- 
tally and theoretically. Theoretical expression for the hyperfine splitting can be calculated 
in the QED framework in the form of a perturbation theory expansion in a, Za, mg/m^. 
Current theoretical uncertainty of this expansion is estimated to be about 70-100 Hz, respec- 
tive relative error does not exceed 2.3 x 10~^ (see discussions in [l-3|). The experimental 
error of the best measurements H, 5| of the muonium HPS is in the interval 16-51 Hz. A 

n 

new higher accuracy measurement of muonium HFS is now planned at J- PARC, Japan [6|. 
Combining muonium HFS theory and experiment one can determine the value of a'^{m^/me) 
with the uncertainty that is dominated by 2.3 x 10^^ relative uncertainty of the HFS theory 
. This is currently the best way to determine the precise value of the electron-muon mass 
ratio. Further reduction of the uncertainty of this mass ratio requires improvement of the 
HFS splitting theory. Main sources of the theoretical uncertainty are due to still unknown 
three-loop purely radiative contributions, three-loop radiative-recoil contribution, and non- 
logarithmic recoil contributions (see detailed discussion in {2, 3|). We consider reduction of 
the theoretical error of HFS splitting in muonium to about 10 Hz as the current goal of the 
HFS theory. 

As a step in this direction we calculate below a three-loop contribution to HFS generated 
by the light by light scattering diagrams in Figs. [1] and [2] with insertions of one loop 
polarization in the upper and lower photon lines, respectively. 
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FIG. 1. 



II. CALCULATIONS 

We start with the light by light scattering contribution to HFS that was calculated long 
time ago It is generated by the diagrams in Fig. [3l where we have not shown explicitly 
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FIG. 2. 

three more diagrams with the crossed photon hnes. In our calculations below we will follow 
the general approach developed in [?! and start with the light by light scattering contribution 
in Fig. |3](see fl) 
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(1) 
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where k'^ is the four-momentum carried by the upper photon lines, q'^ = (0, q) is the spacelike 
four-momentum carried by the lower photon lines, S'^^^" is the light by light scattering 
tensor, and all momenta are measured in the electron mass units. The Fermi energy is 
defined as 



-{Za) (1 + a^) — — 



m, 



(2) 



where is the muon anomalous magnetic moment. The angle brackets in Eq. ([T]) denote 
the projection of the 7-matrix structures on the HFS interval (difference between the states 
with the total spin one and zero). 
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FIG. 3. 



With account for three more diagrams with crossed photon lines not shown explicitly in 
Fig. [3] contributions to HFS of the first two diagrams coincide and we can represent the 
light by light block as a sum of two contributions corresponding to the first two (ladder) 
diagrams in Fig. |3] and corresponding to the crossed (last) diagram in Fig. [3] 



r ^4 

J ir'^i ' 
where (we return to dimensionful momenta here) 
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Calculating traces we obtain 



(4) 



(5) 



X [{p ■q) + {p-k) + {k- q)] + 32^^^" [(A; ■ g)pV - (P ■ qWp'^ - (p ■ k)p''q^] (6) 
+16(p ■k){p- q)g^"g''^ - 32p^p^k^q^} ^ ^\ ^ , 

and 

^a/3Mv _ |8Di^'^" {-•ik^'q^ + fcV + P''/] - SDa^'^^^V - ^Dj,g^''p''q^ 

+16(fc • q)g^''p''p^ - 16(A; ■ q)g>'''k^p^ - 16{k ■ q)g^''p''q^ + 16p ■ (A; + q)g^''k''q^} (7) 
1 

where 



= p'^—m^, D2 = {p—qY—rn?, = {p—kY—m^, D4 = {p—q—kY—m^. (8) 

After calculation of the integrals in Eq. (|3]) we obtain the light by light scattering block 
in the form 

where 
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1) 



1) 



X 



X 



+ 



e{l -y) + q^yz + {k ■ q){2 - y + yz) 
3 2 



+ 



{k\\ - + qS''z^\\ \{k ■ q)g^'^g^^ - 2g^'^k^q^ 



'2{^-y)yz 

Q?{l,y,z, 1) 



rg^ + {k ■ qY 



fi2(l,y,2;, 1) 



(10) 



x{l — x) 

1, 2;, 1) 



Jo Jo 



x{l + x) ^ x{l — x) 



9liaku(ll3 — 
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X y 



n{x,y,l,l) n{x,y,0,l) 
2x'^y 



n{x,y,z, 1 



-{k ■ q)g^,agup + 



9?{x,y,z, 1) 



X 



k^{l - xy) + g^(l - a; + xyz) + (A; ■ g)(l - xy){l - x + xyz) g^aKqp 



(11) 



Vt{x^ y, z, ^) = — k'^xy{l — xy)—q'^x{l — yz){l — x + xyz)—2k-qxy[l—x — z{l — xy)]^, (12) 

and (?^- = (1,-1, -1,-1). 

The 7" matrix structures in Eq. ([T]) are antisymmetric in (a, /3) and (yU, i^), and we have 
thrown away all symmetric in (a, (3) and (/x, i/) terms in Eq. (JTOl) and Eq. (iTTl) . We have also 
combined terms that coincide after antisymmetrization, and deleted even in k and q terms 
that disappear anyway after substitution in the odd in these momenta integral in Eq. ([1]). 
As a result we automatically subtracted symmetric in k and q logarithmically divergent 
contribution in S"''^'^", and the result in Eq. ([9]) is finite and gauge invariant. 

Next we substitute the light by light scattering tensor in Eq. ([T]), and introduce two 
new Feynman parameters t and u to combine the upper photon propagators, the electron 
propagator, and the denominator Q{x, y,z,C,) in the integral representations in Eq. ffTOl) and 
Eq. (ITT]) of the light by light scattering tensor 

(1 - u) [(1 - t)k^ + t{P - 2mko)] + u 



^{x,y,z,^) 
-xy{l - xy) 



{k-Qf-A, 



(13) 
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where Q = grf + r, A = g{—q^ + «^)5 t = ""^(l — u)t, and 



d = 



1 — X 
1 — xy 



u{l-yz){l-x + xyz) _ . 2 _ 1 
1/(1 - 5- 



im?u 



xy{l - xy) 



(14) 

After the Wick rotation and integration over k and q (we return to dimensionless momenta 
here) we obtain an expression for the hght by hght contribution to HFS 



AE = 2AEl + AEc = {2AeL + Aec) ^-Ep, 



where 



(15) 



^eL(C) = EM1c)- 



The integrals Ae^*|^j arise in calculations of the ladder and crossed diagram contributions 
and have the general form 



(16) 



L(C) 



The integrands ^7^*^ and JT^*'' are collected in Tables [T] and [TTl respectively. Notice that not 
all Feynman parameters arise in all integrands in Tables HI and HTl Some parameters just do 
not arise in particular integrals, or take a fixed value, for example, a; = 1 in the ladder light 
by light scattering diagram diagram, see Eq. (fTOj) . As a result the expressions for A in the 
Tables are simpler than the general expression below Eq. f|T3|) . 

The third columns in Tables |T] and [III contain separate integrals Ae^'' and Ae^'' , and 
respective sums in the last lines. The sum of the ladder and crossed diagram contributions 
in Eq. (fTSll nicely reproduces the old result [?, 8] 
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(17) 



AE = -0.472 514 (1 



E, 



(18) 



for light by light scattering contribution to HFS. 

Let us calculate contributions to HFS generated by the diagrams with polarization inser- 
tions in Figs. [1] and_|2l We use the well known integral representation for the polarization 
operator (see, e.g.. 



TABLE I. First Set of Ladder Light by Light Integral^ 



1 

2 
3 

4 
5 
6 

7 



^yz{l-t){l-uf[^ + ^^) 



\z=l 



4 

TT^ 1 — 1/ 



{i-m-u) 



2 / J_ _ r2 
A 



15=1 



4 {l—u)d 
1^ A 15=1 

4 2 (l-t)(l-M)^rf 

i-j, a^ 15=1 



^^(i-t)(i--)^(2k+#;i^^, 
-^^^(i-o(i--)^(2^+^)",^, 

-^(i--Mi-^)(2^ + ^ 



2/3 



2A- 



+ 



2(g2d2_^2) 
A^3 



If) 8y-^ «(1— n)ii 

-•-^ "IF^^^^i; A^ 15=1 
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¥^(l-tMl-n)^d(^ + ^)^^^^ 

12 -¥^a-^Mi-^)^<j^-^),,, 



IC=i 



Ae 
-0.6255. 

0.1498. 
3.2252. 

0.0697. 
0.1628. 
2.0905. 

-3.8178. 

-0.3303. 

-0.4842. 

-0.0193. 
-0.0105. 

-0.0058. 



-0.914. . . 

0.215. . . 
6.403. . . 

0.085. . . 
0.364. . . 
3.927. . . 

-8.070. . . 

-0.669. . . 

-1.078. . . 

-0.036. . . 
-0.019. . . 

-0.010. . . 



0.4045. 



0.195. 



In this table 



d = ^uz, T = (1 — g = 



uz{l-yz) 2 2 
— ] d , a 

1 - y 



y(i - y) 



where the dimensionless momentum q is Euchdean. 

Momentum q in the integrands in Tables [T] and HIl is also Euclidean and to account for 
the polarization operator insertions in both lower photon lines in Fig. [T] it is sufficient to 
insert the factor 2n(g^) in the integrands in Tables HI and UTI Similarly to Eq. f lT5|) respective 
contributions to HPS can be written as 

AEd = 2A£:f + ^E^J = (2Ae7 + Ae^^) ^^^Ep, (20) 



TABLE II. First Set of Crossed Light by Light Integral^; 



1 
2 
3 
4 
5 
6 
7 



J, 



C 



y{l-xy) (-^ *)(-'^ 



(l-t)(l-n)2(i- 

T^yi^-m-u? 



t2 



4 -u) 





y{l-xy) 
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l~x+xyz 




y{l-xyy^ 


4 
_2 


l~x+xyz 



ll/=l 



' A A^ 



I A^ 



^gf^(i-tMi-.)^d( 



A^ 



Ae 



C 



-L6733. .. 
-0.2729. . . 
-0.3665. . . 

-0.2997. . . 

0.3460. . . 

0.9869. . . 
-0.0020. . . 



Ae 



vpii) 
C 



-3.294. . . 
-0.405. . . 
-0.818. . . 
-0.43L . . 

0.387. . . 

2.943. . . 
-0.004. . . 



-1.2816... 



-1.623. . . 



In this table 



d = u 



l-x 
1 - xy^ 



T = (1 — u)t, g 



xy{l ~ xy) 



u{l — yz){l ~ X + xyz) 

2/(1 - y) 

A = g(q2 + a2). 



where 

A4V) = E^^I(C)- (21) 

i 

The fourth columns in Tables HI and HTl contain the integrals Ae^^''*'' and Ae^*'*\ and their 
sums in the last row. Collecting these contributions we obtain the total contribution to HFS 
generated by the diagrams in Fig. [1] 

AErf = -1.2326(5) ^^^f^Ei.. (22) 

Calculation of the contributions to HFS generated by the diagrams in Fig. [2] follows the 
same general route as for the diagrams with polarization insertions in the lower photons. We 
again parameterize contribution to HFS generated the light by light scattering diagrams in 
Fig. [3] exactly like in Eq. ( !T5|) . However, now it is convenient first to integrate analytically 
over momentum q in the integrals in Eq. ([T7|) . As a result the separate contributions to HFS 
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acquire the form 



dy dz dt du d^JC^^y 
'0 Jo Jo Jo Jo 

The integrands JC^l^ and JC^Jj' are collected in Tables [Till and HVl 



TABLE III. Second Set of Ladder Light by Light Integralf|! 



IC 



Ae 



Ae 



vp{i) 



(23) 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 



0-9 ag^ J 



ag ~^ ag'^ 
2 {l-u)d 

15=1 

1 yz (l-t){l-u)^d 



lf=l 



u 



ag op 



2 I 3_ , eP 
ag ag^ 



-^{1- z)u(l-u) i^ + A 

7r V / \ / \ an ' aa^ 



a^g'^ 



3 I 2,if 
op' op" 



2 yz(l—z) u{l—u)d 



7r 1-y 



«5 15=1 



■f^(l-t)n(l-n)2d(^ 
■i^(l-t)n(l-n)2d( 



+ 



d2 
Op 



2 



?=1 



-0.6255. . . 

0.1498. . . 
3.2252. .. 

0.0697. . . 

0.1628. .. 

2.0905. . . 
-3.8178. . . 
-0.3303. . . 
-0.4842. . . 
-0.0193. . . 
-0.0105. . . 
-0.0058. . . 



-0.4631. 
0.1044. 
2.9251. 

0.1083. 

0.1273. 

1.6319. 
-3.3474. 
-0.5174. 
-0.4946. 
-0.0307. 
-0.0168. 
-0.0092. 



0.4045. 



0.0177. 



In this table 



d = £,uz, T = (1 — u)t, g = 



uzjl - yz) _ 2 2 _l 

1 - y 9 1 



2/(1 - y) 



Insertion of the polarization operator in Eq. f lT^ (with — )■ —k^) in the upper photon 
lines of the diagrams in Fig. [2] is described by insertion inside the integrands in Tables IIIII 
and HV] of the factor 
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TABLE IV. Second Set of Crossed Light by Light Integral^; 



1 

2 
3 
4 
5 
6 
7 



C 



1 l+x Q 
27r 



1 



u 



_2_ 



27r y{l-xy) 



r2 



2 = 1 



2=0 



ag 



aFg^ 



1 X 



TT 1—xy 

1 m(I--u) /" _2_ 
7rj;(l-xy) \a9 

1 l—x+xyz/-^ 
2iT y{l—xy)-^ ^ 

1 l—x-\-xyz_/-y 
2-K y{l-xy) y 



ll/=l 



op 



t)u{l — U 
t)u{l - ufd 



2 

0-9 



Ae 



(0 
c 



-1.6733. . . 
-0.2729. . . 
-0.3665. . . 

-0.2997. . . 

0.3460. . . 

0.9869. . . 
-0.0020. . . 



Ae 



vp{i) 
C 



-1.513L .. 
-0.2646. . . 
-0.3084. . . 

-0.2088. . . 

0.5995. . . 

0.8396. . . 
-0.0038. . . 



-1.2816. .. 



-0.8597. . . 



In this table 



l-x 
1 - xy^ 



(1 - %i)t, g 



u(l — yz){l — X -\- xyz) 

y(i ~ y) 



i2 2 

a , a 



xy{l - xy)_ 



introduction of an additional Feynman parameter w and the substitution 



a= aHw) = a- + (25) 

The extra Feynman parameter w does not arise (it is effectively equal one) in the vacuum 
polarization integrals in the 4th, 8th, and 10th rows in Table llllt and in the 5th row in Table 
HVl For these integrals the substitution reduces to — )■ a^(l). 

Similarly to Eq. fl20|) we represent the contributions to HFS generated by the diagrams 
in Fig. |2]in the form 

AEu = 2AEf + AE^J = (2AeF + A^c ) ^^^^f, (26) 

where 

AeI^c) = E^^I(§- (27) 

i 

The fourth columns in Tables UTTl and [TVl contain the integrals Ae^^^*'* and Ae^^*'*, and their 
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sums in the last line. Collecting these contributions we obtain the total contribution to HFS 
generated by the diagrams in Fig. |2] 



AEu = -0.8242(1) "'^^^-' ef. (28) 

71 



III. CONCLUSIONS 

Collecting results in Eq. fl22|) and Eq. fl28l) we obtain the total contribution to HFS 
generated by the polarization insertions in Figs. [T] and H] 



AE,p = -2.056(1)^^^^, (29) 



or numerically 



AE^p = -2.63 Hz. (30) 



Three-loop contribution to HFS containing closed electron loops and factorized one-loop 
radiative insertions in the electron line were calculated earlier jo-lll. Combining those 
corrections with the result in Eq. ( l29l) we obtain the sum of all gauge invariant three-loop 
radiative corrections to HFS calculated thus far 



AEt = -3.338(1) \ ' Ef, (31) 

or numerically 

AEt = -4.28 Hz. (32) 

Work on calculation of the remaining three-loop contributions to HFS is now is progress. 
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